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ATOMIC SPECTRAL FEATURES DURING THERMONUCLEAR FLASHES ON NEUTRON STARS 
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ABSTRACT 

The gravitational redshift measured by Cottam, Paerels and Mendez for the neutron star (NS) in the low-mass 
X-ray binary EXO 0748-676 depends on the identification of an absorption line during a type I burst as the Ha line 
from hydrogenic Fe. We show that Fe is present above the photosphere as long as M > 4 x lO~ 13 M0yr"' during the 
burst. In this limit, the total Fe column is A/p e ~ 3 x 10 19 cm" 2 for incident material of solar abundances and only 
depends on the nuclear physics of the proton spallation. The Fe destruction creates many heavy elements with 
Z < 26 which may imprint photo-ionization edges on the NS spectra during a radius expansion event or in a burst 
cooling tail. Detecting these features in concert with those from Fe would confirm a redshift measurement. We 
also begin to address the radiative transfer problem, and find that a concentrated Fe layer with kT = 1 .2 - 1 .4 ke V 
and column A/p e = 7-20 x 10 20 cm" 2 (depending on the line depth) above the hotter continuum photosphere is 
required to create the Ha line of the observed strength. This estimate must be refined by considerations of non- 
LTE effects as well as resonant line transport. Until these are carried out, we cannot say whether the Fe column 
from accretion and spallation is in conflict with the observations. We also show that hydrogenic Fe might remain 
in the photosphere due to radiative levitation from the high burst flux. 

Subject headings: diffusion - nuclear reactions - stars: abundances, surface - stars: neutron - X-rays: binaries, 
bursts 



1. INTRODUCTION 

The dominant thermal flux during type I bursts and super- 
bursts offer an ideal opportunity for studying neutron star (NS) 
properties through spectral measurements of their photospheres 
(e.g. London, Taam & Howard 1984, 1986; Foster, Ross & 
Fabian 1987). Though never confirmed, the reports by Waki 
et al. (1984), Turner & Breedon (1984), Nakamura, Inoue & 
Tanaka (1988) and Magnier et al. (1989) of 4 keV absorption 
lines during type I bursts motivated theoretical work on the for- 
mation of the Fe Lya line and edge (e.g. Foster et al. 1987; Day, 
Fabian & Ross 1992). Cottam, Paerels & Mendez (2002; here- 
after CPM) observed the accreting NS EXO 0748-676 using 
the Reflection Grating Spectrometer (RGS) aboard the XMM- 
Newton satellite and claimed detection of hydrogen-like and 
helium-like Fe absorption features from the n = 2 — > 3 Ha tran- 
sitions during type I bursts. This Letter is an initial study on the 
implications of the hydrogen-like feature and predictions for 
heavy element abundances near the photosphere. The helium- 
like feature will be left for future work. 

XMM/Newton observed EXO 0748-676 for 335 ks during 
which 28 type I X-ray bursts were seen. CPM combined the 
burst spectra and split them into a high luminosity (L) part 
(from the initial rise and fall) and a low L part, which covers the 
fading to the baseline L. In the high L part, the detected absorp- 
tion line was identified as the Ha transition of hydrogenic Fe. 
In the low L spectrum, a line identified as the He-like Ha tran- 
sition was seen. For these identifications, CPM measured the 
NS's gravitational redshift at 1 +z = 1 .35, giving R = AAGM/c 2 
from the line centroid (neglecting the rotational corrections of 
Fujimoto 1985; Ozel & Psaltis 2003). These spectra have 
photospheric temperatures &r e ff i00 pa 1.6 keV and kT^f^ < 
1 .4 keV, using the measured color temperatures (kTc =1.8 keV 
and kTc < 1 -5 ke V; CPM) and correcting for the redshift and 



spectral hardening from Comptonization (London et al. 1986; 
Madej 1997). 

In § 2, we calculate the Fe abundance when accretion is ac- 
tive, showing that proton spallation of Fe (Bildsten, Salpeter 
& Wasserman 1992; hereafter BSW) naturally yields Np e = 
3.4 x 10 19 cm" 2 for accretion of solar metallicity material. The 
spallation also produces large columns of elements with Z < 26, 
which could be observable via their photoionization edges dur- 
ing the cooling phase of the burst or a radius expansion burst. 
In § 3, we show that Stark broadening (Paerels 1997) domi- 
nates the Ha line equivalent width and that the measured equiv- 
alent width (properly redshifted to the NS surface) of Wh q = 10 
eV implies an Fe column of A/p e ~ 7-20 x 10 20 cm" 2 depend- 
ing on line depth. However, these abundance calculations are 
highly uncertain due to non-LTE (NLTE) effects and resonant 
line transport which we believe play an important role in the 
line formation. We summarize in §4 and note that the radiative 
flux during the burst can levitate hydrogenic Fe. 

2. HEAVY ELEMENT ABUNDANCES FROM ACCRETION 

The continuum photosphere for the 1 .276 keV photons of the 
Fe Ha transition is at an electron column density of N e w avfj, 
where (Tjh is the Thomson cross section. If Fe is uniformly 
abundant at the solar value (/|«3x 10~ 5 by number; Grevesse 
& Sauval 1999) above the continuum photosphere at 1 keV then 
N Fe ks fi/a-rh ~ 5 x 10 19 cm" 2 . Joss (1977) showed that con- 
vection during the burst cannot bring burning fluid elements to 
the photosphere as the entropy there exceeds that in the burning 
layers by large factors. Hence, we look to Fe in the accreted ma- 
terial as the likely source of the photospheric elements. BSW 
studied the deposition and separation of accreting elements, and 
we use their results here. 

Imagine first that accretion halts during the burst. In a 
pure H atmosphere, the upward pointing atmospheric elec- 
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trie field is m p g/2, so that the force on Fe is F = —43m p g. 
The counter-balancing drag force on the drifting Fe nucleus 
is Fdr = Z 2 r ep (m p /m e ) l l 2 pv/nip, where r ep is the electron- 
proton coupling from BSW. The downward drift speed is v « 

1 cm s _1 r 7 3//2 //9 so that it takes t s ~ 1 s(y/l g cm~ 2 )/T^ 2 for 
Fe to fall to a column depth y. Hence, if accretion halts, Fe 
would sink beneath the photosphere in much less than the burst 
duration. 

Given the sub-Eddington peak luminosity, accretion likely 
remains active during the bursts. Kuulkers et al. (2002) have 
shown this to be the case for GX 17+2. The persistent flux 
from EXO 0748-676 when type I bursts occur implies M > 2 x 
10" 10 M Q yr" 1 (Gottwald et al. 1986; CPM). In order to obtain 
a lower limit on the accretion rate per unit area we assume that 
the accretion covers the whole NS, so that m > 10 3 g cm" 2 s" 1 . 
The measured redshift implies R = AAGM/c 2 (CPM), so that 
the NS's surface is inside the last stable orbit regardless of its 
rotation rate (Kluzniak & Wagoner 1985; Kluzniak, Michelson 
& Wagoner 1990). The resulting accretion gap allows for the 
accreting beam to strike the atmosphere with a velocity com- 
parable to free-fall (Kluzniak & Wilson 1991), realizing the 
Coulomb stopping scenario of Zel'dovich and Shakura (1969). 
The integration of a free trajectory from the last stable orbit to 
the surface yields an impact angle sin 9 w 0. 1 and kinetic energy 
of « 200 - 300 Me V/nucleon. 

Since some angular momentum and energy can be lost in the 
accretion gap, we presume that sin 8 > 0.1 and fix the incident 
energy at 200 Me V/nucleon, resulting in stopping a particle 
of charge Ze and mass Am p after traversing a vertical column 
y s i=a 5.4sin6'(mp/(TThM/Z 2 . BSW showed the repercussions of 
the A/Z 2 scaling: (i) the incident electrons stop first, (ii) hy- 
drogen and helium penetrate furthest to y s p , and (iii) all heavier 
elements stop prior to the H/He beam. For example, Fe stops at 
Js.Fe = 0.08y s . p and is spalled by the energetic protons (of cur- 
rent J p w m I nip) until it has sunk to safety (i.e. y > y SjP ). This 
is the mechanism shown by BSW to destroy accreted heavy el- 
ements. Grazing incidence reduces the distance that must be 
fallen to reach safety and increases the surviving fraction of the 
incident CNO elements relative to that in BSW. 

The infalling Fe could have bound electrons as it spirals into 
the NS surface. However, the combination of varying Doppler 
shifts and the short flight time compared to the residence time 
in the atmosphere (derived below) means that the inflow is an 
unlikely environment for narrow line formation. We thus focus 
our discussion on the thermalized Fe in the atmosphere. At 200 
MeV per nucleon, « 30% of the incident Fe suffer a charge 
changing nuclear reaction before traversing « 1 g cm" 2 and 
thermalizing. Accretion is thus depositing thermal Fe into the 
NS atmosphere aty s .F e = 1.0 sin 6* g cm" 2 at arate^Vpe ~ 2/j-/ p /3. 
For the spallation cross section (i.e. that cross section for 
56 Fe to be chipped to lower Z) cr D w 600 mb (Silberberg et 
al. 1998) the Fe lifetime to spallation from the proton beam is 
fd = (vdJ p Y x ~ 2.8 ms(10 3 g cm" 2 s _1 /m). The Fe must sink to 
a depth of y SjP = 14 sin 9 g cm" 2 for safety, so that when m > 

2 g cm" 2 s" 1 (0.1/sin6»)r 7 3/2 (or M > 4 x 1O" 13 M yr" 1 ^ 2 ), 
destruction dominates settling (i.e. t s > fd). 

Hence, as long as M during the burst exceeds 1 % of that oc- 
curring prior to the burst, the destruction by the proton beam 
completely determines the Fe residence time in the upper at- 
mosphere. This seems a condition easily satisfied for many 
bursters and greatly simplifies the calculation of Np e as well 



as the Stark broadening discussed in § 3, since, in this limit, 
the Fe is present only near the stopping point. Because of the 
complicated interplay of the currents and forces in the upper 
atmosphere (see BSW), the local abundance near the stopping 
point ranges from m 1—40 times the incident abundance. The 
proton spallation reduces the Fe abundance as it sinks, imprint- 
ing an abundance profile oc exp(-f s /fd), which is an exponential 
decay with depth. The Fe rich region above the Fe poor layer 
puts a small inverted gradient in the mean molecular weight 
(/i) that would trigger a rapid overturn in an isothermal layer if 
\dlnfi/dlnP\ > 2/5. Gradients this large will not be reached 
as long as the accreting material has an abundance less than ten 
times solar. Doubly diffusive instabilites might still occur, but 
those timescales are likely no shorter than f s . 

Hence, for finding Np e , we only need to balance the Fe de- 
struction rate (neglecting the a beam is less than a 20% effect, 
as those cross sections are rarely larger than twice those induced 
by protons; Ferrando et al. 1988) Np e = Np e aoJ P , with the Fe 
deposition rate N Fe = 2fiJ p /3, giving 

Due to neutron knockout reactions, about 10% of those Fe nu- 
clei are 55 Fe and 5% are 54 Fe. The Fe column is comparable to 
the uniform solar abundance calculation because <jd ~ OTh- 

Since many spallation reactions occur as the nuclei sink, we 
also found the abundances of all lower Z and A nuclei. We cal- 
culated this nuclear cascade by following the evolution of each 
isotope under the constraint that production equals destruction. 
The cross-sections from Silberberg et al. (1998) were used. We 
found that the sum of all nuclei with 20 < Z < 26 outnumber Fe 
by a factor of 2.5-3.0 (see Table 1). In the absence of the BSW 
spallation scenario, these nuclei would not be present. The sin- 
gle most abundant heavy nucleus is 56 Fe. 

The immediate question that comes to mind is whether these 
spallation products can be detected. They may be detected 
via their Ha lines, but the calculation of their predicted EW 
is left to future work. In Table 1, we find the kTin at which 
1/2 of these nuclei would be in the hydrogenic state (for y = 
0.5 g cm" 2 in Table 1), allowing for a possible photoioniza- 
tion edge detection once kT < kT^n (presuming LTE). To fit 
all Z, we used the explicit cross-section from Rybicki & Light- 
man (1979) and set the bound-free Gaunt factor to g = 0.8, 
giving <7 bf (E) = 6.32 x 10" 18 cm 2 (£ , e /£ , ) 3 Z" 2 where £ e is the 
edge energy (shown in Table 1, value in parentheses is the 
redshifted edge energy) from Bethe & Salpeter (1957; here- 
after BS). This Obf agreed to within (3-5)% to Verner et al.'s 
(1996) analytic fits for Ar, Ca and Fe. For now, we estimate 
the impact of such an edge on the spectrum by considering 
the integral over the optical depth above the photoionization 
edge, r e = Nubf, assuming an effectively cold atmosphere (i.e. 
no reemission above the edge) EW e = J(l -exp(-N<Jbf(E)))dE . 
For r e « 1, the EW? integral is simply expanded to obtain 
EW e = 175 eV(AT/10 19 cm" 2 ) integrated from E = E e to 1.3E e 
which is independent of Z. 

Table 1 shows EW e assuming the element is present at the 
level implied by the spallation scenario and finding the amount 
of it in the Hydrogenic ground state from Saha equilibrium for 
kT = 1.2 keV. As the burst cools or during a radius expan- 
sion burst, these edges could become prominent, and the entry 
of Te imax in Table 1 presumes all of the spallation created ele- 
ments have at least one bound electron. It is intriguing to note 
that these edge energies and EW e are close to those previously 
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reported by Waki et. al. (1984), Turner & Breedon (1984), 
Nakamura et al. (1988), Magnier et al. (1989) and Kuulkers et 
al. (2002). 

3. ESTIMATED IRON ABUNDANCES FROM OBSERVED 
EQUIVALENT WIDTHS 

We now estimate the Fe column needed to give the Ha line 
observed by CPM. Motivated by our accretion scenario, we ini- 
tially presume that this line forms in a thin isothermal layer 
of Fe in LTE at temperature T above the hotter continuum 
photosphere. The observed redshift implies a surface grav- 
ity of g « 3 x 10 14 cm s" 2 (1.4M /M). Though rapid rotation 
(s=s 100 Hz) might explain the observed line width (rs 40 eV), the 
line's equivalent width depends most strongly on the intrinsic 
broadening mechanisms. The scale for thermal Doppler broad- 
ening is AE D /E = ^kT/ArripC 2 « 1.5 x 10" 4 (/tr/l keV) 1 / 2 . 
The natural line width, r « 5 x 10 13 s" 1 (BS), is given by the 
spontaneous decay rate for the 3 — > 1,2 transitions, which in 
units of the thermal Doppler shift is a = HT/2AEd ~ 0.3. 

We now show that Stark broadening is larger by an order 
of magnitude than either of these effects for the Ha transi- 
tion. Stark broadening results from the shift of line energies 
due to the electric field from neighboring ions, e/r^, where ro = 
(47rn c /3)- 1 / 3 is the mean ion spacing (Mihalas 1978). The en- 
ergy perturbation from the linear Stark effect is AE$ = e 1 Axjr\, 
where Ax is the size of the orbital. Using the linear Stark effect 
is justified since its energy shift at the densities and tempera- 
tures considered exceeds that of the degeneracy breaking Lamb 
shift for the relativistic H-like ion (BS). For the n = 3 state of 
Fe at a Thomson depth unity, we find AEs/Eo = 1.5 x 10~ 3 , an 
order of magnitude larger than thermal Doppler. 

Since the majority of the intrinsic line width is attributed to 
Stark broadening, N Fe must be found from the measured EW 
given by the curve of growth (Mihalas 1978), 

f°° ri(a) 

W = 2A AE T J TT-Ja, (2) 
Jo 1 + »?(«) 

where W is the EW in energy units, Ao is the normalized line 
depth at line center, AE is the intrinsic energy width (in our 
case the Stark shift), r\ is the ratio of the line opacity to the 
continuum opacity and a = E/AE. We use the electric field 
distribution for a uniform ion sea (e.g. the Holtsmark dis- 
tribution, Wh', Mihalas 1978), giving r](a) = rjoWnia) where 
Vo = ^VFe,n=2fiine/A r cO'Th is the ratio of line optical depth to con- 
tinuum optical depth at line center. Depending on density either 
thermal Doppler or Stark broadening dominates, so we con- 
volve the Holtsmark distribution with the Voigt profile to cover 
the complete range of densities. We also ignore Debye screen- 
ing effects on the electric field distribution since the deviation 
from the Holtsmark distribution is small for the typical T and 
P- 

Because rotational broadening acts to lower the apparent 
line depth and is likely to be important, the observed value of 
Aq « 0.3 is a lower limit. Indeed, the true line depth could be 
significantly larger than 0.3. However, such a large value of Aq 
is difficult to understand for realistic temperature gradients and 
presuming LTE. For a gray atmosphere (a reasonable approxi- 
mation for the temperature profile of a type I burst atmosphere) 
Ao attains a maximal value of 0.2 when the Fe resides at very 
low column depths. The resolution of this dilemma may be the 
role of NLTE effects and resonant line transport in the Ha line. 
Resonant line transport could generate such an absorption fea- 
ture with a large line depth without invoking large temperature 



contrasts as in LTE (Jefferies 1968). NLTE effects clearly play 
a role in determining abundances from the EW measurements. 
For instance, an immediate measure of the role of NLTE ef- 
fects is that the radiative ionization rate for the Fe ground state 
exceeds the collisional ionization rate (London et al. 1986), 
complicating the ionization balance for Fe. 

To resolve the role of NLTE effects, the radiative transfer 
equation must be solved self-consistently for both the resonant 
transport of the Ha and Lya lines with detailed balanced of the 
Fe ionization and excitation in the strong radiation field of a 
Type I burst. For now, we simply present results for Ao « 0.3, 
which is set by observations and Ao = 0.2, which is characteris- 
tic of a gray atmosphere. 

We calculate the Afce required to reach Wn a = 10 eV presum- 
ing all Fe resides at the column depth, y = P/g. We assume 
LTE and that the dominant states for Fe are either fully ionized, 
or in the hydrogenic ground state or helium-like ground state 
to approximate the partition function. For radial infall, with 
the Fe layer between kT = 1.2- 1.4 keV and A = 0.3 (0.2), 
A^Fe = 7 (20) x 10 20 cm" 2 . This temperature range is consistent 
with the absence of helium-like Fe in the high L phase of the 
burst (CPM). Since the occupation of the n = 1 state exceeds 
that in n = 2 by ~ 10 2-3 , we also calculate W^ya for the im- 
plied A^Fe values. For this line, thermal Doppler and natural line 
broadening dominate Stark broadening (Madej 1989). When 
kT = 1.2—1.4 keV and the Fe resides at the pressure implied 
by radial accretion, W^ya < 30 eV. Foster et al. (1987) also 
calculated Wi,y a for a uniformly distributed Fe abundance at the 
solar value and found that it would not exceed 30 eV. It would 
require an Ao ~ 1 to give an Fe column consistent with that 
calculated from spallation. However, such large temperature 
contrasts (and hence such large values of Ao) are not possible 
in these atmospheres and we need to understand the NLTE line 
transport problem before we can address this apparent incon- 
sistency. 

Day, Fabian & Ross (1992) pointed out that a uniform pres- 
ence of Fe at the solar value would create an observable photo- 
ionization edge. The Fe column require to produce the observed 
Ha line would yield enormous edges. However, NLTE effects 
would weaken the strength of these photoionization features 
significantly (London et al. 1986). The reported broad spec- 
tral features in the 6-9 keV range (van Paradijs et al. 1990, 
Strohmayer & Brown 2002) during type I bursts and super- 
bursts need to be re-analysed as possible Fe photo-ionization 
edges. 

4. SUMMARY AND CONCLUSIONS 

Active accretion in the context of a beam hitting the atmo- 
sphere naturally gives Afce ~ 3 x 10 19 cm -2 which only depends 
on the physics of proton spallation of Fe. In this scenario, heavy 
elements with Z < 26 are produced in quantities comparable 
to Fe and may be observed via their photoionization edges or 
atomic lines (see Table 1) during the cooling phase of type I 
bursts or superbursts. Our initial LTE calculations have found 
that an Fe column of w 7-20 x 10 20 cm" 2 is needed to produce 
the measured Wu a of CPM. However, until a more complete 
calculation of the resonant line transport for the Ha and Lya 
line has been undertaken in this strong radiation environment, 
we cannot say whether this factor of twenty discrepency is cer- 
tain. 

Throughout this paper, we have presumed that Fe is present 
because of active accretion during the burst. However, the ra- 
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diative force on hydrogenic Fe in the ground state is F ra( j = 
/ (JeFe/c dE, where F E is the flux per photon energy at the 
Lya transition (E = 6.9 keV). Since F E is nearly constant across 
the line, F m & = ^he 2 /i^Fe / 'm e c 2 in the optically thin limit 
(Michaud 1970), where /i_>2 = 0.42 is the oscillator strength 
from BS. If F E is given by a black-body, then at kT e s = 1 keV, 
the radiative acceleration on Fe is a rdd « 10 15 cm s" 2 , a factor of 
three larger than g. The radiative force will be even larger for 
realistic type I burst spectra, which are harder than a blackbody 
at the observed T eB (London et al. 1986; Pavlov, et al. 1991). 
The Fe diffusion time is less than the burst duration, so there 
is adequate time for the accelerated Fe to reach a new diffu- 
sive equilibrium. The calculation of the resulting Fe column, 
which could well saturate the line and lower the radiative force, 
is beyond the scope of this Letter, as it requires a simultaneous 



solution of the radiative transfer and the Fe diffusion equations, 
like that done for hot white dwarfs (e.g. Chayer, Fontaine & 
Wesemael 1995). 
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Table 1 



Heavy Element Abundances and Atomic Features 





E e 


7~e,max 


Ni/Ne, 


EW e 






[keV(redsh)] 






[eV] 


[keV] 


Fe 


9.279(6.87) 


0.32 


1.00 


277 


1.12 


Mn 


8.573(6.35) 


0.18 


0.36 


66 


1.05 


Cr 


7.897(5.85) 


0.26 


0.48 


56 


0.99 


V 


7.248(5.37) 


0.26 


0.43 


32 


0.93 


Ti 


6.627(4.91) 


0.31 


0.48 


22 


0.86 


Sc 


6.035(4.47) 


0.41 


0.58 


17 


0.81 


Ca 


5.471(4.05) 


0.47 


0.59 


11 


0.75 


K 


4.934(3.65) 


0.48 


0.55 


7 


0.69 


Ar 


4.427(3.28) 


0.53 


0.54 


4 


0.65 



